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Assemblies of magnetic nanoparticles exhibit interesting physical properties arising from the com- 
petition of intraparticle dynamics and interparticle interactions. In ordered arrays of magnetic 
nanoparticles magnetostatic interparticle interactions introduce collective dynamics acting competi- 
tively to random anisotropy. Basic understanding, characterization and control of dipolar interaction 
effects in arrays of magnetic nanoparticles is an issue of central importance. To this end, numer- 
ical simulation techniques offer an indispensable tool. We report on Monte Carlo studies of the 
magnetic hysteresis and spin-dependent transport in thin films formed by ordered arrays of mag- 
netic nanoparticles. Emphasis is given to the modifications of the single-particle behavior due to 
interparticle dipolar interactions as these arise in quantities of experimental interest, such as, the 
magnetization, the susceptibility and the magnetoresistance. We investigate the role of the struc- 
tural parameters of an array (interparticle separation, number of stacked monolayers) and the role 
of the internal structure of the nanoparticles (single phase, core-shell). Dipolar interactions are re- 
sponsible for anisotropic magnetic behavior between the in-plane and out-of-plane directions of the 
sample, which is reflected on the investigated magnetic properties (magnetization, transverse suscep- 
tibility and magnetoresistance) and the parameters of the array (remanent magnetization, coercive 
field, and blocking temperature). Our numerical results are compared to existing measurements on 
self-assembled arrays of Fe-based and Co nanoparticles is made. 
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I. INTRODUCTION 

Laterally confined magnetic nanostructures (dots and 
nanoparticles) is an important class of novel materials 
with unique physical properties, that emerge because 
their size becomes comparable to various characteristic 
physical lengths (correlation length, domain wall width, 
etc). Owing to their novel physical properties they find 
numerous technological applications in magnetic storage 
media 1 , magnetic sensors^ and magnetic logic devices^ 

Magnetic nanoparticles (NPs) are commonly formed 
in assemblies, with either random or ordered structure. 
In the first group belong systems such as ferrofluids and 
granular solids, while in the second group belong the pat- 
terned media (or magnetic dots) and the self-assembled 
arrays (SAA) of NPs. The existence of order in a NP 
assembly is a decisive property in view of their applica- 
tion in magnetic storage media with ultrahigh density 
(~ Tb/in 2 ), which rely on the possibility of treating the 
NPs as individually addressable magnetic bits. 

Magnetic dots are developed by lithographic processes 
and they are characterized by lateral dimensions in the 
range of 10 — 100 nmA- The major advantage of this ap- 
proach to development of ordered nanostructured mate- 
rials is that it offers great flexibility and good control over 
the shape, the size and the arrangement of the dots as 
well as the choice of the constituent material. However, 
shape imperfections and formation of polycrystalline dots 
are the major factors that determine the deviations from 
perfect periodicity in the arrays. Furthermore, the size of 
the magnetic dots, determined by the limitations of the 



lithographic process is in most cases comparable to the 
exchange length, thus permitting the domain formation 
in the ground state or during magnetization reversal, a 
fact that makes the study of the magnetic behavior of 
dot arrays quite intricate. Growth methods and mag- 
netic properties of isolated dots and dot arrays have been 
recently reviewed by Martin et alA 

In magnetic NPs with diameter D « 1 — 10 nm, the 
domain wall width is well beyond the diameter of the 
NP, and consequently domain formation at the ground 
state or during magnetization reversal is precluded. The 
magnetization of the NPs is practically at its saturation 
value up to temperatures less than but close to the bulk 
Curie temperature. They are commonly referred to as 
single-domain (SD) particles, to distinguish them from 
magnetic dots that exhibit inhomogeneous magnetiza- 
tion. Highly monodisperse (a < 5%) SD magnetic NPs 
are produced by solution chemistry methods and orga- 
nized in hexagonal ordered arrays by self-assemblyi^'^ 
The small size, the high monodispersity and the array 
periodicity achieved by self-assembly, as well as the low 
production cost of SAA, has motivated a great deal of re- 
search effort in the field of synthesis and magnetic charac- 
terization of these systems. The synthetic routes and the 
structural and magnetic characterization methods have 
been recently reviewed by Willard et al£. 

Assemblies of magnetic nanoparticles have been in- 
vestigated in an effort to gain basic understanding of 
the interplay between single-particle magnetic anisotropy 
and interparticle magnetostatic interactions. To this 
end, a variety of sample preparation methods have 
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been adopted by research groups in order to grow 
and magnetically characterize nanoparticle assemblies, 
as for example, frozen ferrofluid o 10 ' 11 ' 12 , discontinuous 
metal-insulator multilayer s 13 ' 14 ' 15 ' 16 , co-sputtered metal- 
insulator film a ' 18 ' 19 , cluster-assembled film o 20 ' 21 self- 
organized particle arrays on surfaces^ 2 - and chemically 
produced self-assembled nanoparticle arraysi 5 ' 6 ' 23 ' 24 ' 25 ' 26 
Investigations of the static and dynamic magnetic prop- 
erties of dipolar interacting nanoparticle assemblies 
brought up fundamental issues related to the existence 
of a ground state which shares common features with 
canonical spin glasses (slow relaxation, memory and age- 
ing effects) i 12 ' 13 i 14 ' 15 i 27 . Dipolar interparticle interac- 
tions are considered responsible for the observed com- 
plex (spin-glass-like) behavior of sufficiently dense and 
random nanoparticle assemblies. In view of the tech- 
nological importance of the hysteresis behavior and the 
thermal stability of magnetization of NP assemblies 2 ^, 
the important issue of the effects of magnetostatic inter- 
actions on the static magnetic properties has also been 
studied extensivel y 10 ' 11 

Contrary to random assemblies, ordered arrays of NPs 
are ideal systems to investigate the role of interparticle 
interactions, for two reasons. First, the NP arrange- 
ment is periodic with small perturbations, and the unde- 
sired complications introduced by spatial randomness are 
substantially suppressed. Second, chemically synthesized 
magnetic NPs are often coated by an inorganic surfactant 
layer that prevents agglomeration during self-assembling 
but also keeps the surfaces of neighboring particles at a 
distance well beyond the range of exchange forces. As a 
consequence, the prevailing interparticle interactions in 
a SAA are magnetostatic. Finally, the spherical, in most 
cases, shape of the NPs diminishes the importance of 
higher order multipolar interactions and the assembly is 
well described by dipolar interparticle forces. 

Various experiments have demonstrated the presence 
of magnetostatic interactions in SAA with various de- 
grees of structural disorder and layered NP assem- 
blies. Reduction of the remanence at low temperature,— 
increase of the blocking temperature ; 18 ' 30 ' 31 increase 
of the barrier distribution widthj 3 ^ deviations of the 
zero-field cooled magnetization curves from the Curie 
behavior,^ 3 - difference between the in-plane and out- 
of-plane remanence,— and increase of the blocking tem- 
perature with frequency of applied field 3 ^ have been ob- 
served and attributed to interparticle dipole-dipole inter- 
actions (DDI). Long-range ferromagnetic order in linear 
chain o 22 ' 36 ' 37 and hexagonal array s 24 ' 25 ' 26 ' 34 of dipolar 
coupled single-domain magnetic nanoparticles has been 
demonstrated, supporting the existence of a dipolar su- 
perferromagnetic ground state. Nanoparticle assemblies 
with random morphology have been studied more as most 
growth techniques developed so far (sputtering, clus- 
ter beams, mechanical alloying) produce random sam- 
ples. Ordered nanoparticle arrays, on the other hand, 
have been less studied both theoretically and experimen- 
tally due to the difficulty in producing ordered samples. 



Chemical synthesis and self-assembly offer a new and 
promising approach to this direction.— We therefore be- 
lieve that basic understanding of the magnetic properties 
of dense (interacting) ordered arrays is currently highly 
demanded. 

In the ongoing research effort for development of mag- 
netic nanostructures with reduced size and improved 
thermal stability^ 8 -, the exploitation of the exchange bias 
effect in laterally confined structures (dots and nanopar- 
ticles) has attracted a lot of interest^ Atomic scale mod- 
els of the magnetic structure have been developed in an 
effort to interpret experimental observations of the ex- 
change bias effect in composite NPs with a ferromag- 
netic (FM) core and an antiferromagnetic (AFM) shell. 
Among the most important theoretical results^ - we men- 
tion (i) the disappearance of the exchange bias field (He) 
at temperatures above the Neel temperature of the AFM, 
in agreement with experiments^ 1 -, (ii) the strong depen- 
dence of He on the number unsaturated bonds across 
the FM-AFM interface and the dependence of He on 
the interface area, (iii) the increase of both He and He 
for a given NP radius with increasing oxidation depth, 
(iv) the increase of He and decrease of He with increas- 
ing oxidation layer thickness and a fixed core radius, (v) 
the fast stabilization of He with increasing core size, in 
agreement with experiments^ 2 -, and (vi) the reduction of 
He and increase of He and its thermal stability with in- 
creasing exchange constant of the AFM material and/or 
at the FM-AFM interface. Despite the research effort 
focused on the atomic scale mechanism of magnetization 
reversal in composite nanoparticles , 40 ' 43 ' 44 ' 45 much less 
attention has been paid so far to the modification of the 
magnetic hysteresis behavior due to inter-particle inter- 
actions arising in assemblies. In this direction, Fe NPs 
embedded in iron-oxide matrix— were shown to freeze 
below a temperature owing to the competition between 
the exchange anisotropy at the core-shell interface and 
the interparticle DDI. Similarly, increase of the exchange 
bias field due to magnetostatic interparticle coupling was 
found in stripes of Co/CoO nanoparticles^ 7 - and inter-dot 
magnetostatic interactions were shown to produce asym- 
metric anomalies in the magnetization reversal mecha- 
nism of Co/CoO dot arrays^ 8 - The modification of the 
coercive and exchange-bias fields in dense nanoparticle 
arrays with core-shell morphology as a result of the com- 
petition between exchange anisotropy and interparticle 
dipolar interactions consists a challenging issue. 

Detection and quantification of DDI in assemblies 
of magnetic NPs has been addressed so far by a 
variety of experimental techniques including in most 
cases SQUID magnetometry and AC susceptibility 
measurements^ -, and more recently small-angle neu- 
tron scattering (SANS)— and resonant magnetic X-ray 
scattering^ -. The last method is a direct probes of mag- 
netic correlations at the interparticle scale. These studies 
have provided ample evidence that the interplay between 
random anisotropy and DDI determine the magnetic be- 
havior of the NP assemblies. More recently dipolar in- 
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teraction effects and the resulting collective dynamics in 
SAA of FeSi, Fe 2 0^, and C<M^- NPs was studied by 
reversible transverse susceptibility (RTS) measurements. 
The RTS technique is a well established and powerful 
method to obtain information about the anisotropy of 
magnetic nanoparticles, from considerations of the peak 
positions of the field-dependent RTS^ The implemen- 
tation of RTS to study SAA arrays revealed information 
regarding the different dynamical regimes of an interact- 
ing assembly accessed as the temperature increases. 

Electron spin is a degree of freedom whose control 
and detection in transport measurements is the basis 
of the rapidly developing field of spintronics^ 6 -. Charge 
transport measurements in a SAA of Co NPs were 
performed 5 -^ and revealed a spin-dependent tunneling 
mechanism which is responsible for substantial (~ 10%) 
tunneling magnetoresistance (TMR) values at low tem- 
perature (~ 20K). The tunneling barriers are provided 
in an array by the insulating surfactant layer surround- 
ing the NPs. TMR measurements probe the interparticle 
correlations within the range of the spin-diffusion length, 
and are therefore sensitive to the magnetic microstruc- 
ture of the assembly. Interparticle interaction effects are 
expected to reveal themselves in the TMR signal. The 
field-dependent magnetization and conductivity were dis- 
cussed in the experiments of Black et al^, however a sys- 
tematic correlation between the two quantities remains 
to be performed. 

In addition to the experimental work, various numer- 
ical studies that focused on the ground state configu- 
ration and the hysteresis behavior of dipolar interact- 
ing nanoparticle arrays have appeared. The interplay of 
DDI and perpendicular anisotropy was shown^ to in- 
duce a reorientation transition below a critical temper- 
ature and interaction-induced shape anisotropy of a fi- 
nite sample controls the magnetization reversal mode. 
Dipolar interactions were found to decrease the coer- 
cive field of magnetic nanoparticle arrays independently 
of the array topology (square or hexagonal) despite the 
fact that the ground state configuration is determined by 
the array topology^ The presence of an incomplete sec- 
ond layer with hexagonal structure does not destroy the 
long-range FM ordering of the ground stated while even 
slight structural disorder within the array destroys that 
ordering^ On the other hand, higher order (quadropo- 
lar) magnetostatic interactions were shown to act in syn- 
ergy with DDI stabilizing the long range order of the 
ground state in a nanoparticle array^ Previous theoret- 
ical studies of RTS in random assemblies of magnetic NPs 
demonstrated that a wide size distribution rounds the 
peaks of RTS 63 , orientational texture suppresses the co- 
ercivity peak^ 3 - and dipolar interactions lead to merge of 
the coercivity and anisotropy peak s 64 ' 65 . More recently 
the issue of the structure of the RTS curves of SAA was 
addressed by Monte Carlo (MC) simulations^ 6 - that re- 
produced many of the experimental observations from 
RTS measurements in SAA of Fe NPs^i. Charge trans- 
port in nanoparticle arrays has been studied by resistor 



network models (RN) that include in a phenomenological 
way the essential aspects of the thermally-activated hop- 
ping mechanism 6 ^, the spin-dependence of the hopping 
proccess 6 ^ and the details of the micromagnetic configu- 
ration of the sampled In a recent study 6 ^ the signature 
of dipolar interaction effects in TMR measurements has 
been investigated. 

In this article we review our results from MC simu- 
lations of the field and temperature dependence of the 
magnetization, the RTS and the TMR of hexagonal ar- 
rays of dipolar coupled magnetic nanoparticles with ran- 
dom anisotropy. The consideration of a hexagonal ar- 
rangement of NP is an essential feature of our model as 
DDI have a well known anisotropic character that relates 
their magnitude and sign to the relative position of the 
interacting dipoles. The main structural parameters we 
focus on are (a) the interparticle separation, which is di- 
rectly related to the dipolar coupling strength and can 
be experimentally controlled by variation of the surfac- 
tant layer during the synthetic process, and (b) the sam- 
ple thickness, namely the number of stacked monolayers 
(MLs) , which is a crucial parameter for the collective re- 
sponse of the array, controlled by the NP concentration 
in the colloidal dispersion. The aim of the present work is 
to reveal the modification introduced to experimentally 
measured properties of an ordered NP assembly (mag- 
netization, susceptibility, magnetoresistance) due to the 
presence of DDI. The remaining of the paper is orga- 
nized as follows: In Section II we present the structural 
and magnetic model used in our simulations. In Section 
III we discuss numerical results on the hysteresis charac- 
teristics (saturation remanence, coercivity), the zero field 
cooled (ZFC) and field cooled (FC) curves and and the 
extracted blocking temperature (Tb). The magnetization 
of the interesting class of composite nanoparticles with a 
core-shell morphology is also discussed. The evolution of 
the RTS curves with temperature and dipolar strength 
is discusses next, and finally magnetoresistance calcula- 
tions are presented. Whenever experimental results are 
available they are compared with our simulations aiming 
to reveal the character of interparticle interactions in dif- 
ferent measured samples. Final conclusions and remarks 
are given in the Section IV. 



II. MODEL AND SIMULATION METHOD 

We proceed with the definition of the spin model used 
to describe the magnetic structure of a magnetic NP ar- 
ray formed by either simple ferromagnetic or composite 
(FM core/ AFM shell) nanoparticles. The NPs form- 
ing an ordered array are assumed spherical and monodis- 
perse. The size dispersion is not expected to introduce 
major modifications to the magnetic behavior because 
particularly low values are achieved (a < 5%) in most 
samples! 6 ' 9 ' 25 The NP diameter is D and they occupy 
the sites of a triangular lattice in the zy-plane with lat- 
tice constant d > D. When more that one MLs are con- 
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sidered, the particles are close-packed in an ABC ABC... 
stacking sequence. This structure is consistent with elec- 
tron microscopy studies of Co & 30 ' 33 and FeQ NP arrays. 
Incomplete layer of NPs are formed by random occu- 
pancy of the triangular lattice sites. The NPs are single 
domain with uniaxial anisotropy in a random direction, 
and they interact via dipolar forces. The total energy of 
the system is 



\ Si • Sj 3{Si ■ Rij)(Si ■ Rij) 

E = g l^ Rj. 

ij »j 

-fc£($.$) a - A (1) 



where Si is the magnetic moment direction (spin) of the 
i-th particle, el is the easy-axis direction, and Rij is the 
center-to-center distance between particles i and j. Hats 
in Eq (1), and further on, indicate unit vectors. The en- 
ergy parameters entering Eq. |T]) are the dipolar energy 
g = m? I 'd 3 , where m — M 8 V is the particle moment, 
the anisotropy energy k = K{V, and the Zeeman en- 
ergy h = mH due to the applied dc field H. The energy 
parameters (g, k, h) entering Eq. |T]), the thermal energy 
t = ksT, and the history of the sample determine the mi- 
cromagnetic configuration at a certain temperature and 
bias field. Because, our simulation method relies on min- 
imization of the free energy of the system, multiplication 
of all the energy parameters by the same scaling fac- 
tor does not modify the results. Thus, in all subsequent 
results we scale the energy parameters entering Eq. ([T]) 
by the single particle anisotropy energy (k = 1). This 
choice makes our numerical results applicable to a class 
of materials with the same parameter ratios rather than 
to a specific material. The crucial parameter that deter- 
mines the transition from single-particle to collective be- 
havior is the ratio of the dipolar to the anisotropy energy 
g/k = {■K/Q)(M 2 JK 1 )(D/df. The reported values^i^ 
for fee or hep Co NPs are g/k = 0.2 - 0.4(D/d) 3 , while 
for the soft e-Co phase, higher values are expected^ For 
Fe NPs Farrell et a£Z. report g/k = 1.54 {D/df. Despite 
the relative dispersion of the reported values, the impor- 
tant issue is that for most samples of Co and Fe NPs the 
ratio of dipolar to anisotropy strengths is below unity, 
except for e-Co. Thus, in the numerical results presented 
in the following section we consider g/k values less than 
one. 

Extensions to the spin model described by Eq. ([I]) are 
required in order to study composite NPs with a FM 
core and an AFM shell. In the present work we adopt a 
model introduced by Meiklejohn and Bean 7 ! (further on 
referred to as the MB model) in their interpretation of 
shifted loops observed in oxidized transition metal NPs 
samples after zero-field cooling. The MB model provides 
a phenomenological understanding of the exchange bias 
effect and the unidirectional anisotropy! 7 -^ Consequently, 
important parameters of the exchange-bias effect such as 
the interface structure and interface magnetization are 
averaged out. Despite its simplicity, the MB model and 



its variations was successfully implemented in the case of 
FM/AFM bilayers to interpret the dependence of the ex- 
change bias field on temperature^ 3 -, on the thickness the 
AFM layer— and on the direction of the applied field. 75 . 
The major weaknesses of the MB model being the over- 
estimation of the exchange-bias field values^ and the 
underestimation of the coercivity values^ We adopt the 
MB model as the simplest possible approach to bring out 
the essential aspects of the competition between intra- 
particle (uniaxial anisotropy) and interparticle (dipolar) 
interactions. Our purpose, is to investigate this inter- 
play, which is expected to be important in dense samples, 
rather than revealing the atomic scale mechanism which 
is responsible for the exchange bias effect. 

According to the MB model, coherent rotation of the 
atomic spins is assumed in the FM core and the AFM 
shell, while the net magnetic moment of the shell is van- 
ishingly small. In addition, the interface of the AFM is 
assumed fully uncompensated, namely all spins belong to 
the same sublattice and is exchanged coupled to the core. 
Consequently, the magnetic state of each NP is described 
by a pair of anisotropic and exchange coupled spins, S FM 
and S AFM . For a dipolar interacting assembly the total 
energy reads 

qFM qFM tfeFM 3 \(qFM 3 ^ 

E = 9 2^ rJ- 

-kc £(Sf M ■ S) a ~ k s J2(Sf FM ■ Sif 

i i 

Sf M ■ Sf FAI -hJ2sf M -H (2) 



where Sf M and Sf FM indicate the magnetization direc- 
tions (spins) of the core and the interface layer of the 
shell, respectively. J is the interface exchange energy 
and kc , ks are the distinct values of the core and shell 
anisotropy, respectively. Notice that owing to the zero 
net magnetization of the AFM shell, only the FM cores 
couple to the external field and between them via magne- 
tostatic forces. For simplicity we have assumed in Eq. (0 
that the core and shell magnetizations have a common 
easy axis, that is therefore labeled by the particle in- 
dex i. Due to the vanishing net magnetization of the 
shell there is neither Zeeman nor dipolar contributions 
to the total energy due to the shell. However, the cou- 
pling across the FM-AFM interface makes an exchange 
contribution to the total energy, expressed by the fourth 
term in Eq. @. 

For either simple or composite NP assemblies, the mag- 
netic configuration is obtained by a MC simulation, us- 
ing the standard Metropolis algorithm^ The initial spin 
configuration corresponds to the saturation state along 
a chosen axis (x or z). Experimentally, to observe a 
shifted loop a field-cooling process is performed prior to 
the hysteresis measurement in order to align the AFM 
moments parallel to the moments of the FM. 39 i 72 Fur- 
thermore, the value of the observed exchange bias field 
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FIG. 1: Sketch of the sample and applied magnetic field used 
in our simulations, (a) In-plane and (b) out-of-plane direc- 
tions of the bias field (Hdc)- The indicated (H ac ) field is con- 
sidered only in the calculation of the transverse susceptibility. 
For the magnetoresistance calculations a dc bias voltage is ap- 
plied on opposite edges of the sample along the j/-axis. 



the array. For a given micromagnetic configuration {Si} 
we introduce the spin-dependent conductivity between 
two nanoparticles i and j as£&, 

(Jij = itq(1 + P 2 cos9ij) exp(—Rij/a - E c /k B T) (4) 

where <jq = 2e 2 /h is the conductivity quantum, P is the 
spin polarization of the conduction electrons, cosOij = 
(Si ■ Sj) , E c — e 2 /2C is the the activation energy to 
charge a neutral NP by addition of a single electron, C is 
the NP capacitance relative to its surrounding medium, 
and a = hj \/Sm*(U — Ef) is the electron wave function 
decay length in the insulating barrier of height U relative 
to the Fermi energy. In all our simulations we take a = 
d, as a sufficient requirement to allow charge transfer 
between neighboring nanoparticles and P = 0.34 which 
is an appropriate for Co NPs i 30 ' 83 Charge conservation 
on every node of the network implies 



i.e. the loop shift) increases with the value of the cooling 
field (gfc) ' 39 ' 72 The choice of the saturation state as the 
initial state to calculate the hysteresis loop of core-shell 
NPs from the fully saturated state, it is equivalent to as- 
suming an infinitely strong cooling field (Hfc ^ Jafm)- 
Thus, the maximum value of He is obtained. During re- 
laxation, the initial 10 3 MC steps per spin (MCS) are 
used for relaxation of the system towards equilibrium 
and thermal averages are calculated over the subsequent 
10 4 MCS, allowing 10 MCS between sampling events to 
achieve statistical independence. The results are aver- 
aged over N c = 30 — 100 samples with different realiza- 
tions of the random axes distribution and the thermal 
fluctuations. To deal with the long-range character of 
the DDI we use periodic boundaries in the xy-plane and 
implemented the Ewald summation method adapted to 
a quasi-two-dimensional system / 8 ' 79 ! 80 Free boundaries 
along the z-axis are assumed. 

Transverse susceptibility measurements are performed 
with an ac measuring field (H ac ) perpendicular to the dc 
bias field (see, Fig. [1]). The ac field is weak (~ 10 Oe) 
and its frequency lies in the rf regime (/ ~ 10 6 Hz) i 52 ' 81 . 
The weak measuring field permits us to neglect transverse 
hysteresis effects and calculate the RTS in the zero-field 
limit. Furthermore, since that the Neel relaxation time 
of the NP magnetization is large compared to the in- 
verse frequency of the ac field the static approximation 
for the measuring field is justified. Susceptibility values 
are obtained from the fluctuations of the magnetization 

xl {±) (H x{z) ) = -j^[{M 2 y ) (My) 2 } (3) 

where x^r^ is the in-plane (out-of-plane) RTS and N is 
the number of NPs in the simulation cell. 

In the last part of this section we describe the resistor 
network model used to calculate the magnetoresistance of 



^°ij{<t>i - fa) = ® (5) 

ij 

where fa is the local value of the electric potential. 
Eq. ([5]) is solved for the local potentials with the bound- 
ary conditions that set the local potential on opposite 
sides of the sample, namely at y — and y = L, (see 
Fig. Q] equal to zero and fa, respectively. The total con- 
ductivity is given as 

ij 

Obviously the values of a depend on the spin configu- 
ration Si. A thermal average of the conductivity is ob- 
tained by averaging the conductivity values, as obtained 
from Eq. ^ over a sequence of equilibrium spin config- 
urations produced by the MC algorithm. Finally, the 
tunneling magnetoresistance of the sample is defined as 

TMR(H) = Gsat ~ a N (g) (7) 
a(H) 

where a sa t is the saturation value of the conductivity. It 
follows from Eq. Q that the local conductivity between 
particles i and j increases quadratically with spin po- 
larization (P) and exponentially with localization length 
(a) and particle capacitance (C), the latter depending 
on the NP diameter.— The TMR values are expected to 
have a similar dependence on the parameters a and P. 
However, for a monodisperse sample the TMR values are 
independent of C. Since we are mainly interested here 
is the shape of the field-dependent TMR curves, rather 
than the actual values of TMR, the dependence of the 
TMR on a and P is not considered further on. 
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Dipolar Strength (g) 

FIG. 2: Dependence of (a) saturation remanence and (b) coer- 
civity on the dipolar strength at low temperature (t/k = 0.01) 
for a monolayer of FM nanoparticles. Closed symbols : In- 
plane field. Open symbols : Out-of-plane field 



III. NUMERICAL RESULTS AND DISCUSSION 

A. Isothermal hysteresis and ZFC/FC 
magnetization of FM nanoparticles 

Two important characteristics of the hysteresis loop 
are the saturation remanence (M r ) and the coercivity 
(H c ). The former is the magnetization of the system af- 
ter removal of a saturating field, and the latter the nega- 
tive field required to zero the magnetization of a system 
being in the positive remanent state before application 
of the field. In Fig. [2] we show the dependence of M r and 
H c on the strength of the dipole-dipole interaction. Vari- 
ation of the dipolar strength in a SAA can be achieved 
experimentally in two different ways, namely, by varying 
the thickness of the surfactant layer—, or by synthesis of 
different material NPs with the same size and the same 
surfactant layer thickness, or combination of both. In 
Fig. [5] we show that depending on the direction of the ap- 



plied field these quantities show different variation with 
dipolar strength. This behavior can be explained by the 
ferromagnetic and anisotropic character of the DDI. It is 
a well establishe d 59 ' 62 ! 85 fact that DDI on a triangular 
lattice stabilize a FM ground state and create an easy- 
plane for the array magnetization due to their anisotropic 
character. Therefore, the in-plane remanence tends to- 
wards the saturation value (M r /M s = 1) as the dipolar 
strength increases, while the out-of-plane remanence de- 
creases continuously because the out-of-plane field is nor- 
mal to the easy-plane. On the other hand, the coercivity 
decreases with increasing dipolar strength, independently 
of the applied field direction. This behavior can be un- 
derstood as due to collective rotation of the moments. 
For an in-plane field, the effective anisotropy of a cluster 
of dipolar coupled NPs is reduced (as a result of an av- 
eraging process over many random easy directions) and 
the total moment of a cluster is larger than a single NP. 
Due to the synergy of these two factors a weaker reversal 
field is required. For an out-of-plane field, the develop- 
ment of an easy-plane forces the moments to lie in the 
xy-plane reducing their projection along the field axis, 
thus a weak field is required to fully zero the magnetiza- 
tion. Anisotropy between the in-plane and out-of-plane 
remanence (Mj > M^~) has been observed in arrays of 
Co NPs 34 , where a ratio of 7 = M^/mJ = 0.30 was 
found. For these samples the reported^ 4 - dipolar strength 
is g/k — 0.07 and from the data shown in Fig. [2] we 
obtain -f(g = 0.07) = 0.3 in good agreement with the 
experiments. A zero-temperature calculation 34 on a tri- 
angular lattice gave a similar value for 7. In the same 
experiments^ 4 - a negligibly small dependence of the coer- 
civity on the applied field direction was found. Our sim- 
ulations agree with this feature showing, that for dipolar 
strengths in the range g/k = — 0.1, Hq and Hq practi- 
cally coincide (Fig.[T]). Our MC simulations and the en- 
ergy minimization approach 34 justify the dominant role 
of DDI in the magnetic properties of these arrays. 

An important structural parameter in SAA is 
the number of stacked MLs. This parameter can 
be partially controlled either by varying the sol- 
vent concentratio n 8 i 23 i 30 i 33 prior to self-assembly or 
by forming the arrays using the Langmuir-Blodgett 
tcchniquei 35 ! 86 To address the effects of film thickness on 
the hysteresis properties of NP arrays we show in Fig. [3] 
the variation of the remanence and coercivity with the 
number of MLs. For strong dipolar coupling the rema- 
nence takes maximum values when complete monolayers 
form, while an incomplete top layer suppresses the mag- 
netic order due to the competing character of the DDI 
in a random nanoparticle assembly 87 . Peaks in H c are 
observed at full coverage, which indicate that fully or- 
dered samples are magnetically harder relative to sam- 
ples with disordered surfaces. The observed decay of H c 
with increasing layer thickness marks a transition from a 
two-dimensional reversal mode (c = 1ML), during which 
the moments are forced by the interaction-induced easy- 
plane anisotropy to remain in-plane during reversal, to a 
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FIG. 3: Dependence of saturation remanence (upper panel) 
and coercivity (lower panel) on number of stacked monolay- 
ers of FM nanoparticles, at low temperature (t/k = 0.01). 
Non-integer values of film thickness correspond to a ran- 
domly occupied uppermost layer. Triangles: weak coupling 
(g/k = 0.25). Circles: strong coupling (g/k = 10). The ap- 
plied field lies in-plane. 



three-dimensional mode (c > 2 ML), during which the ro- 
tation path of the moments is not restricted in the plane 
of the film. When DDI are weak, it is shown in Fig. [3] 
that the oscillatory dependence of M r is suppressed and 
the peaks in H c are washed out. However, the decay 
of H c with increasing thickness remains, reaching a con- 
stant value above c ~ 2 MLs. The increase of M r with 
coverage, observed for submonolayer coverage in Fig. [3J 
is in agreement with experiments on dilute samples of 
Fe NPs^S, with variable concentrations. Also SAA of Fc 
NPs always showed higher M r values compared to dilute 
(disordered) samples^. This trend is reproduced by our 
simulation results in Fig. [3J 

As temperature increases thermal fluctuations can as- 
sist the magnetic moments to overcome the anisotropy 
barrier leading the system to thermal equilibrium. The 
regime in which this is achieved is defined by the blocking 



FIG. 4: Temperature dependence of (a) saturation remanence 
and (b) coercivity of a monolayer of weakly dipolar (g — 0.1) 
FM nanoparticles. Closed symbols: In-plane field. Open 
symbols: Out-of-plane field. Dashed line: Non-interacting 
nanoparticles. 



temperature of the system. For dipolar interacting NPs 
the concept of a single-particle barrier becomes rather 
vague as the thermal activation of a moment is correlated 
to the motion of all the other moments of the system. 
However, one can still refer to the blocking temperature 
of the system in a phcnomenological way, namely one can 
define it as the temperature above which the remanence 
and coercivity vanish. We show in Fig. Q] the temperature 
dependence of the in-plane and out-of-plane M r and H c 
for a monolayer of weakly coupled NPs. We first notice 
that the effect of interactions is to increase the blocking 
temperature of the array (tf, ~ 0.17) relative to the non- 
interacting case (tf, ~ 0.14). This is clear in Fig. 0] for the 
in-plane field. Most interestingly, at temperatures above 
t;J, M r and H c of the interacting system are non zero. 

This result defines an interesting temperature regime, 
t£ < t < t;,, in which the thermal energy overwhelms 
the nominal anisotropy barrier (Ef> ~ fc), but the hys- 
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FIG. 5: Temperature dependence of the first nearest-neighbor 
correlation function at the remanent state for a monolayer of 
weakly dipolar (g — 0.1) FM nanoparticles. Closed symbols: 
In-plane field. Open symbols: Out-of-plane field. Dashed 
line: non-interacting assembly. 

teresis behavior of the NP assembly persists due to DDI. 
Calculations of the energy barrier distribution in dipo- 
lar interacting NP assemblies showed that, DDI broaden 
the distribution thus producing a range of high barri- 
ers (Ef, > k), which are responsible for the observed en- 
hanced thermal stability in dense assemblies. 88 Notice 
finally, that the enhanced thermal stability is less pro- 
nounced, but also existent in the out-of-plane geometry, 
despite the fact that in this geometry dipole interactions 
have a demagnetizing effect at low temperature (Fig. [5]) . 
A similar enhancement of M r and H c at elevated tem- 
peratures was previously predicted also for NP assemblies 
with random morphology 87 and its existence was verified 
experimentally from the measured increase of the block- 
ing temperature with NP density in frozen ferrofluids.— 
Further insight into the collective behavior at elevated 
temperatures can be extracted from examination of the 
short-range moment correlation function, defined as 

SRCF = (Si ■ S 3 ) Rtj=d - (St) 2 . (8) 

Strictly speaking, the blocking temperature of a sys- 
tem corresponds to the maximum value of the long range 
moment correlation function, namely the susceptibility. 
We prefer to examine the short-range correlation func- 
tion because it is directly accessible by various spectro- 
scopic experiments to be mentioned below. The position 
of the SRCF peak of gives a reasonable approximation 
to the blocking temperature of the system. The temper- 
ature dependence of the SRCF, shown in Fig. shows 
that short range FM correlations exist above the blocking 
temperature of the non-interacting array (t? ~ 0.14) and 
persist up to temperatures above the blocking tempera- 
ture of the interacting array (tjj ~ 0.17). Recently, Kor- 
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FIG. 6: ZFC/FC magnetization (at h/k — 0.1) for a mono- 
layer of weakly interacting (g = 0.1) FM nanoparticles. 
Closed symbols: In-plane field. Open symbols: Out-of-plane 
field. Dashed lines indicate the different values of the in-plane 
and out-of-plane blocking temperatures. 

tright et a/ 50 extracted the interparticle magnetic corre- 
lations in dense arrays of Co NPs from x-ray scattering 
experiments. They concluded that for e-Co NPs, which 
are strongly dipolar, AFM correlations exist at temper- 
atures above the blocking. SANS studies of Fe NP as- 
semblies showed no evidence of AFM correlations at el- 
evated temperature. 89 Finally, our simulations indicate 
the existence of short-range FM interparticle correlations 
in dipolar coupled arrays. There seems to be difficult at 
present to reach a decisive conclusion about the nature 
of interparticle correlations in NP arrays from experi- 
mental findings. Scattering experiments (x-rays, SANS) 
at lower Q-values, are expected to probe magnetic cor- 
relations at a scale lying well within the inter-particle 
separation regime are most probably required in order to 
compare with our present simulations 

An interesting feature shown in Fig. [5] is the different 
temperatures at which the peaks of SRCF are located. 
The relative shift of the peaks indicates that the array 
of dipolar interacting NPs exhibits anisotropic blocking 
behavior with t} > Tf. 

The anisotropic blocking can also be demonstrated 
by examination of the ZFC/FC curves shown in Fig. [5] 
The existence of anisotropic blocking for dipolar inter- 
acting arrays of Fe NPs has been recently demonstrated 
experimentally in dense arrays of Fe NPs^i In these 
experiments a ratio tJ'/T^ ~ 1.15 was obtained from 
ZFC/FC measurements. The Fe NP parameters given^i 
are M s — 1360 emu/cc, K — 3.4-10 4 erg/cc, D = 6.8 nm 
and d ~ 20 nm, which lead to g/k ~ 0.11. Our simu- 
lations shown in Fig. [5] for g/k — 0.1 give t'i/t^ ~ 1.08, 
which is in reasonable agreement with the experiments, 
given various factors not considered in our model, such as 
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FIG. 7: (a) Dependence of blocking temperature (a) on dipo- 
lar strength, for a monolayer of FM nanoparticles, and (b) on 
film thickness, for FM nanoparticles with g/k = 0.25. Mea- 
suring in-plane field h/k — 0.1. The values of Tj, are normal- 
ized to the value corresponding to non-interacting nanoparti- 
cles. 



the sample thickness, the deviations from perfect stack- 
ing of the monolayers and the in-plane structural defects. 

Changing the structural parameters of a self-assembled 
array, namely interparticle spacing and film thickness is 
expected to affect the blocking behavior. The blocking 
temperature is proportional to the effective barrier height 
which is determined by the single-particle anisotropy and 
the local dipolar field. Dipolar interactions on the other 
hand are sensitive to interparticle spacing (g ~ 1/d?). 
Decreasing the interparticle spacing is equivalent to in- 
creasing the dipolar strength, and consequently modify- 
ing the blocking temperature. Similarly, increasing the 
film thickness we shown above that it suppresses at low- 
temperature the values of M r and H c . In Fig. [7] we show 
the effect the structural changes have on the blocking 
temperature of the assembly, the latter being obtained 
from the peak of the corresponding ZFC curve. We see 
that in the range of dipolar strengths considered here, Tf, 



scales linearly with the dipolar strength, or equivalently, 
it decreases with the cube of the interparticle separation 
(T(, ~ 1/d 3 )- We can therefore think of the effect of weak 
dipole interactions (g/k < 1) as an increase of the single- 
particle anisotropy barrier by an amount proportional to 
their coupling strength. 

A similar effect on 7], to decreasing the interparticle 
separation is obtained by increasing the areal coverage 
and the film thickness. As shown in Fig. [7] a dramatic 
increase of T}, is observed during formation of the first 
monolayer and a saturation behavior is reached as soon 
as the second complete monolayer is formed. Reduction 
of Ti, upon dilution of chemically synthesized assemblies 
has been reported by several group a 6 ' 25 i 30 ' 31 . In partic- 
ular, Zhang et al^- performed ZFC/FC measurements 
on self-assembled e-Co NP and report 30% increase of 
X}, relative to highly dilute samples. The NP parame- 
ters for these samples^ correspond to g/k ~ 0.2, thus 
the predictions of our simulations (Fig. UJ are in agree- 
ment with these experiments. Similar dependence of T& 
on NP concentration has been also observed for NPs dis- 
persed in a solid matrix^ and reproduced by simulations 
for 3D disordered assemblies of magnetic NPs. 87,90 . It is 
interesting that despite the demagnetizing character of 
the dipolar interaction in the ground state of 3D random 
assemblie o 87 ! 91 , they tend to stabilize the FM character 
of the assembly at elevated temperatures (T < T h °). The 
increase of X& with number of stacked monolayers was 
observed in discontinuous C0-AI2O3 multilayers 18 . How- 
ever, the slow saturation of Tf> obtained after 5-7 mono- 
layers in these experiments is probably due to deviations 
from the ideal stacking sequence and the in-plane ran- 
domness in NP size and location, inherent to the sam- 
ple preparation technique. We also believe that the lack 
of perfect stacking is probably the reason that measure- 
ments of T& in Langmuir-Blodgett films of Co NPs showed 
a decrease of Tf, with increasing thickness^ 



B. Isothermal hysteresis of FM/AFM nanoparticles 

We discuss next the hysteresis behavior of a 2D hexag- 
onal array of magnetic NPs with core-shell morphology. 
To calculate the hysteresis loop, the spins are set initially 
in the saturation state along the positive x-axifM- and the 
field is swept from positive to negative values and back. 
This choice produces a negative shift of the hysteresis 
loop, namely a positive exchange bias field. The main 
issue we address is the dependence of the remanence, the 
(effective) coercivity and the exchange bias field on the 
strength of the DDI. The total energy of the NP assembly 
is described by EqO The coercivity for the shifted loops 
is defined as He = (l/2)\Hci — Hc%\, and the exchange 
bias field as H E = (l/2)\H C i+H C2 \, where H C i (H C2 ) is 
the upper (lower) branch coercivity corresponding to the 
backward (forward), with respect to the exchange bias 
field direction, magnetization reversal process. 

In the MB model the microscopic details (atomic struc- 
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FIG. 8: Hysteresis loops at low temperature (t/kc = 0.01) of 
a monolayer of core-shell NPs with different dipolar strengths 
and interface exchange values, (a) J/kc = 1.5 , (b) J/kc ~ 
1.0 and (c) J/kc = 0.5. Squares: g = 0, circles: g/kc = 0.3, 
and triangles: g/kc = 1.0. Shell anisotropy ks/kc = 5.0. 
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FIG. 9: Dependence of coercivity (he) and exchange bias 
field (He) on dipolar strength for a monolayer of core-shell 
NPs at low temperature (t/kc = 0.01). Curves for different 
values of the interface exchange are shown. Triangles: J/kc = 
1.5, circles: J/kc = 1.0, and squares: J/kc = 0.5. Shell 
anisotropy ks/kc = 5.0. 
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FIG. 10: Dependence of the upper branch coercivity (hci) 
and the lower branch coercivity (hc2) on dipolar strength for 
a monolayer of core-shell NPs at low temperature (t/kc = 
0.01). Curves for different values of the interface exchange 
coupling are shown. Triangles: J/kc = 1.5, circles: J/kc = 
1.0, and squares: J/kc = 0.5. Shell anisotropy ks/kc = 5.0. 



ture, defects, exchange coupling strength) are absorbed 
into the value of the exchange constant J. Thus results 
for different values of J are considered. The anisotropy 
of the AFM oxide is assumed much higher than the core 
anisotropy, and in the present work we take ks/kc = 5.0. 
This value of ks is high enough to ensure blocking of the 
shell magnetization (Safm) for applied fields in the range 
that the core magnetization (Sfm) exhibits hysteresis 
behavior. Typical hysteresis loops for dipolar interacting 
NP arrays are shown in Fig. [5] Similarly to the case of 
simple FM NPs discussed in the previous section, M r in 



core-shell NPs increases and He decreases with increas- 
ing DDI strength (Fig. [9]). However the changes of He 
due to dipolar coupling are controlled by the value of the 
interface exchange. In particular, larger reduction of He 
is observed in systems with larger interface exchange J. 
Nevertheless, the fraction of He reduction is nearly the 
same for all values of J. 

The behavior of He with increasing dipolar strength 
appears more complex (Fig. [9]). For strong interface cou- 
pling (J/kc = 1-5) He drops linearly with the dipolar 
strength and for weak interface coupling (J/kc = 0.5) 
it is only weakly dependent on the dipolar strength. 
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The most striking behavior is observed when the in- 
terface exchange is comparable to the core anisotropy 
(J/kc = 1-0). In this case, weak DDI enhance He, which 
exhibits a maximum value around g/k c ~ 0.2. The ob- 
served enhancement of He due to weak DDI is contrary 
to what one expects on intuitive grounds. Namely, DDI 
that lead to symmetric backward and forward magneti- 
zation reversal processes, are shown to enhance the loop 
asymmetry. Insight into the dependence of the exchange 
bias field on dipolar strength can be obtained from sep- 
arate examination of the two branches of the hysteresis 
loop and the corresponding coercivity values Hex and 
Hqi- As shown in Fig. [TUl Hex decreases (in absolute 
value) with dipolar strength, independently of the inter- 
face exchange value. This behavior is similar to what is 
observed in FM NPs (Fig. \2§ at low temperature. The 
similarity stems from the fact that in the case of back- 
ward magnetization reversal the exchange bias field in- 
creases the barrier height for reversal of an isolated mo- 
ment. The collective reversal induced by DDI facilitates 
the reversal process leading to (absolute) lower coerciv- 
ity values. On the contrary, the dependence of Hci on 
dipolar strength varies significantly with the interface ex- 
change. In the forward reversal process the exchange 
bias field reduces the barrier for the reversal of an iso- 
lated moment. When the value of the exchange is weak 
(J/kc = 0.5, Fig. fTU]). a reduced but finite barrier exists 
for the forward reversal. As previously, the collective re- 
versal, induced by DDI, reduces further the barrier height 
leading to smaller Hc2 values. When the interface ex- 
change is strong (J/kc = 1.5, Fig. HO]) , the barrier for 
the forward reversal of an isolated moment disappears 
and the forward reversal becomes a "downhill" process 
in the energy landscape, a fact reflected in the negative 
value of Hc2- In this case, DDI introduce additional bar- 
riers due to their anisotropic character and the collective 
motion of the moments obstructs the reversal. The val- 
ues of Hc2 approach zero (from negative values) as the 
dipolar strength increases, reflecting the increasing dif- 
ficulty for forward reversal. For intermediate interface 
exchange values (J/kc — 1.0, Fig. ITOjl the dependence of 
Hc2 on the dipolar strength shows a varying behavior. 
For isolated moments (g = 0), Hc2 assumes small posi- 
tive values indicating a very low barrier for forward rever- 
sal. DDI initially facilitate the forward reversal, but with 
increasing dipolar strength the barrier vanishes and Hc2 
reaches a negative value (g/kc — 0.2). Further increase 
of the dipolar strength inhibits the forward reversal. The 
observed dependence of He on DDI strength shown in 
Fig. [5] is the net effect from the variation of the upper 
and lower branch coercivities with dipolar strength. Ex- 
perimentally, enhancement of He due to DDI has been 
observed^ in stripes of Co/CoO NPs where the quasi- 
one-dimensional morphology of the array enhances the 
effect of DDI. 
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FIG. 11: Field dependence of reversible transverse suscepti- 
bility of a monolayer of FM nanoparticles at low-temperature 
(t/k = 0.1). The bias field is swept from negative to positive 
values. Closed circles: in-plane field. Open circles: out-of- 
plane field. The arrows in the out-of-plane data indicate the 
position of the coercive field. 



C. Transverse susceptibility of FM nanoparticles 

Measurements of the field dependent RTS have been 
long ago predicted theoretically^ to reveal direct infor- 
mation on the magnetic anisotropy of a NP assembly. 
A typical RTS curve of a NP assembly with random 
anisotropy, obtained at low temperature by sweeping the 
bias (dc) field in one direction, exhibits three character- 
istic peaks, one at the coercive field and two at positions 
corresponding to the anisotropy field (±Hk)- The lat- 
ter is related to the single-particle anisotropy as Hk = 
2K\/M S . Thus measurement of RTS should, in principle, 
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FIG. 12: Dependence of anisotropy field on dipolar strength 
at low temperature (t — 0.05) for a monolayer of FM nanopar- 
ticles. Closed circles: in-plane field. Open circles: out-of- 
plane field. 



be a direct method to measure K\. However, wide parti- 
cle size distributions 6 - 3 -, thermal fluctuation s 64 ' 65 ! 66 and 
most importantly, dipolar interaction effect s 64 ! 65 ' 66 ! 93 
modify the position and shape of the peaks making the 
determination of the anisotropy strength uncertain. We 
discuss here the evolution of the field dependent RTS 
curves under increasing values of the dipolar strength, 
or equivalently decreasing values of the interparticle dis- 
tance in SAA of magnetic NPs. We show in Fig. [11] 

the in-plane (xt) an d out-of-plane (Xt) RTS curves 
at low temperature (t/k = 0.05) for increasing dipolar 
strength. The non-interacting sample shows clearly the 
theoretically predicted three peaks located at the coer- 
cive (h c ) and anisotropy (±hx) fields. Downshifted val- 
ues of h c and hx relative to the zero-temperature values 
(h c = 0.98 and Hk — 2) are due to thermal fluctuations 
effects. The most important effects of dipolar interaction 
on the RTS curves, shown in Fig. [TT] and rj2] are : (i) 
the suppression of the h c peak of Xt > (h) the location of 
the Kk peak of Xt a * higher fields than the correspond- 
ing peaks of x\i (iii) the downshift (upshift) of the Kk 
peak of (Xt) with increasing dipolar strength, (iv) the 
slower saturation with bias field of Xt relative to Xt- in 
recent RTS measurement in Fe NP arrays 5 - 1 -, observations 
similar to our points (ii) and (iv) were made. However, 
the coercivity peak could not be resolved probably due 
to not sufficient lowering of the temperature or due to 
the presence of NP size distribution. Nevertheless, the 
agreement of our results in points (ii,iv) above, consti- 
tute sufficient evidence that DDI were responsible for the 
observed experimental trends in these measurements. A 
physical interpretation of the observed opposite trends of 
the in-plane and out-of-plane Kk peaks with increasing 
dipolar strength (Fig. [T2"|) relies on the development of an 
easy-plane anisotropy induced by DDI. For an in-plane 
bias field, the interaction-induced anisotropy reduces the 



barrier for an irreversible switching of the moments lead- 
ing to a reduction of the anisotropy field, while in the out- 
of-plane geometry DDI inhibit the irreversible switching 
along the z-axis by developing an easy-plane normal to 
this axis, thus increasing the anisotropy field. Finally, the 
linear dependence of the anisotropy peaks on the dipo- 
lar coupling, or equivalently, on the inverse cube of the 
interparticle spacing (Hk ~ l/^ 3 ), shown in Fig. [T2l for 
g/k < 0.4, could be used to perform an extrapolation pro- 
cedure on measurements taken at different interparticle 
separations in order to extract the value of single-particle 
anisotropy (K\). 



D. Tunneling magnetoresistance in FM 
nanoparticle arrays 

Tunneling magnetoresistance refers to a large decrease 
of a sample's resistivity under application of a bias 
magnetic field, observed when charge carriers transmit 
through two FM regions separated by a non-magnetic 
insulating barrier. The basic mechanism underlying the 
effect is spin-dependent scattering of the carriers. The 
first FM region acts as a polarizer for the electron spin 
while the second region causes scattering whose strength 
is proportional to the misalignment of the magnetization 
relative to the first region. In the case considered here, 
the magnetic NPs are the relevant FM regions and the 
surfactant layer separating them is the insulating bar- 
rier. Thus, in principle, spin-dependent transport mea- 
surements should reflect the underlying micromagnetic 
structure of the NP assembly. In Fig. [T3] we plot the 
field dependent TMR of a monolayer of dipolar coupled 
NPs and compare it with the corresponding branch of the 
magnetization hysteresis loop. The sharp peak of TMR 
occurs very close to the coercive field, because the spin 
disorder in the array is maximized at his field. The effects 
of DDI can be observed in the TMR curves. A down- 
shift of the TMR peak position with increasing dipolar 
strength is observed following the reduction of the H c 
values. The value of TMR at the remanent state de- 
creases with interactions, reflecting the increasing align- 
ment of the magnetic moments demonstrated also by the 
increased values of the M r (see Fig.f2]). Finally, the TMR 
sensitivity, namely the slope of the field dependent TMR, 
increases with increasing interaction strength, due to a 
collective reversal of the moments during which the de- 
gree of alignment is higher for stronger interactions. 

The easy-pale anisotropy induced by DDI in a 2D-array 
of NPs is expected to produce a strong dependence of the 
TMR values on the direction of the bias field. Indeed the 
TMR curves shown in Fig. [TH vary substantially with the 
azimuth angle (6) . The TMR sensitivity decreases as the 
field approaches the z-axis, reflecting the slow saturation 
of the magnetization for an out-of-plane bias field. A 
striking feature occurring for 9 ~ 15° in Fig. fTJ] is that 
the peak of TMR does not occur at the coercive field 
but at a higher field. We state that the peak of TMR 
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FIG. 13: Field-dependence of the magnetization (upper 
panel) and the corresponding tunneling magnetoresistance 
(lower panel) of a monolayer of FM nanoparticles, at low 
temperature (t/k — 0.01). The curves correspond to different 
dipolar strengths. Circles: g = 0, up-triangles: g/k = 0.1, 
down-triangles: g/k = 0.2, and diamonds: g/k — 1. The 
bias field lies in-plane and is swept from negative to positive 
values. 



FIG. 14: Field-dependence of the magnetization (upper 
panel) and the corresponding tunneling magnetoresistance 
(lower panel) of a monolayer of dipolar interacting (g/k = 0.2) 
FM nanoparticles, at low temperature (t/k = 0.01). The 
curves correspond to different directions of the bias field rel- 
ative to the z-axis (9 = 0°— 90°). The field is swept from 
negative to positive values. 



occurs at the critical field (H Q ), namely the field for an 
irreversible switch of the magnetization^ rather than at 
the coercive field. This has been verified by simulations 
in purely dipolar arrays^ or non interacting arrays with 
aligned easy axes. In both cases the critical field can be 
obtained analytically and it was found that the peak of 
TMR occurs exactly at this field. When DDI are absent, 
the random anisotropy leads to H c « 0.96-ffd2i and the 
TMR peak occurs very close to the coercive field. How- 
ever, as DDI increase, they induce coherent rotation of 
the moments and a dominant easy-plane anisotropy. In 
the case of an easy-plane anisotropy the difference be- 
tween the H a and H c is maximum for 6 « 15°^, which 
explains why the maximum deviation between the field 
corresponding to the peak of TMR and the coercive field 
occurs at this angle. 



IV. FINAL REMARKS 

We have given an overview of our simulation studies 
on the effects of dipolar interactions on the magnetiza- 
tion, the short-range correlation function, the transverse 
susceptibility and the conductivity of nanoparticle as- 
semblies forming hexagonal arrays. We discussed the 
modifications introduced to the field and temperature 
dependence of the above quantities due to the dipolar 
coupling between nanoparticles. We have demonstrated 
that dipolar interactions in thin layers formed by self- 
assembled NPs induce an easy-plane anisotropy that is 
responsible for an anisotropic magnetic behavior with 
respect to the bias field direction. The anisotropic be- 
havior is revealed in different physical quantities used 
to magnetically characterize nanoparticle arrays, such as 
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the transverse susceptibility and the tunneling magne- 
toresistance and in the values of characteristic magnetic 
parameters, such as the remanence, the coercivity and 
the blocking temperature. The most important dipolar 
interaction effects predicted by our numerical studies and 
observed in experiments are summarized as follows : (i) 
The increase of remanence with particle density for sub- 
monolayer coverage, observed in self-assembled Fe NPs^ 
(ii) The enhancement (suppression) of the in-plane (out- 
of-plane) remanence at low temperatures, in agreement 
with magnetization measurements on self-assembled Co 
NPs. 2,5 (iii) The suppression of the coercivity at low tem- 
peratures with decreasing interparticle spacing, observed 
in self-assembled Fe NPsj2£ 

(iv) The increase of the apparent blocking tempera- 
ture with decreasing interparticle separation (Tj, ~ \/d 5 ), 
in agreement with ZFC magnetization measurements on 
self- assembled e-Co NPs^i (v) The existence of two dis- 
tinct blocking temperatures for the in-plane and out- 
of-plane geometries, in agreement with ZFC magnetiza- 
tion and transverse susceptibility measurements on self- 
assembled Fe NPsi^i (vi) The persistence of short-range 
FM correlations at temperatures above the blocking ex- 
istence; x-rays measurements on Co NPs 50 showed AFM 
correlations at high temperatures, while SANS measure- 
ments on Fe NPs did not provide evidence of AFM 
correlations, 4 ? More refined experiments are required in 
this direction, (vii) The possibility to enhance the ex- 
change bias field with decreasing interparticle separa- 
tion at low temperatures in qualitative agreement with 
the enhanced exchange bias field found in oxidized Co 
NPs forming chain-like structures4i (viii) The decrease 
(increase) of the anisotropy field, obtained from trans- 



verse susceptibility measurements with an in-plane (out- 
of-plane) bias field, with decreasing interparticle separa- 
tion, in agreement with measurements on self-assembled 
Fe NP samples £i 

Certain results of our simulations that are yet to be 
verified by experiments include : (i) The transition from 
a 2D to 3D magnetization reversal mechanism with in- 
creasing number of monolayers that is responsible for a 
reduction of the coercivity with number of stacked mono- 
layers. A layer-by-layer measurement of the magnetic 
properties would address this point, (ii) The suppression 
of the TMR values and the increase of TMR sensitivity in 
dense arrays. Verification of this result requires compar- 
ison between TMR measurements in dilute samples and 
self- assembled arrays, and finally (iii) The location of the 
TMR peak not at the coercive field, for a tilted bias field 
in dense assemblies New developments in self-assembly 
techniques^ will make feasible further investigation of 
the interplay between geometry and interaction effects 
in magnetic NP arrays, while multiscale computational 
schemes*^ constitute a promising tool to achieve more 
detailed understanding and control of their properties. 
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